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Abstract. The Sunyaev Zel'dovich (SZ) effect from galaxy clusters is one of the most powerful cosmological tools 
for investigating the large-scale Universe. The big advantage of the SZ effect is its redshift independence, which is 
not the case for visible and X-ray observations. It allows us to directly estimate the cluster's total mass from the 
integrated comptonization parameter Y, even for distant clusters. However, not having a full knowing intra-cluster 
medium (ICM) physics can affect the results. By taking self-similar temperature and density profiles of the ICM 
into account, we studied how different ICM morphologies can affect the cluster total mass estimation. With the 
help of the high percentage of cool core (CC) clusters, as observed so far, the present analysis focuses on studying 
this class of objects. A sample of eight nearby (0.1 < z < 0.5) and high- mass (M > 10^^ Mq) clusters observed by 
Chandra was considered. We simulated SZ observations of these clusters through X-ray derived information and 
analyzed the mock SZ data again with the simplistic assumption of an isothermal 6eia-model profile for the ICM. 
The bias on the recovered cluster total mass using different sets of assumptions is estimated to be 50% higher 
in the case of hydrostatic equilibrium. Possible contributions to the total bias due to the line-of-sight integration 
and the considered ICM template are taken into account. The large biases on total mass recovery firmly support, 
if still necessary, cluster modeling based on more sophisticated universal profiles as derived by X-ray observations 
of local objects and hydrodynamical simulations. 

Key words. Sunyaev Zel'dovich effect - cluster of galaxies - cool cores - individual clusters: Abell 1413, Abell 1689, 
AbeU 1835, Abell 2204, Abell 2261, MS1358.4+6245, RXJ1347.5-1145, ZW3146. 

1. Introduction 

Clusters of galaxies are the largest gravitationally-bound 

. objects arising thus far from the process of hierarchi- 

' cal structure formation (jVoit 20051) . As the most recent 

, and most massive objects of the Universe, clusters are 

■ excellent probes for studying its formation and evolu- 

, tion. The observed state of gas within a cluster is de- 

' termined by a combination of shock heating during ac- 

, cretion, radiative cooling, and thermal feedback produced 

' by the cooling itself, so the density and temperature of 

, the ICM represent the full thermal history of clusters' 
formation. To better understand the physics of ICM, it 
is necessary to have sufficient knowledge of the gas den- 
sity and temperature distributions. Though clusters are 
the ideal target objects for X-ray observations of the hot 
ICM, millimeter and sub-millimeter measurements pro- 
vide independent and complementary tools for studying 
the same ICM by exploiting the Sunyaev Zel'dovich (SZ) 
effect dSunyaev k Zel'dovich 19721 ). 
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The SZ effect is the Comptonization of the cosmic mi- 
crowave background (CMB) photons, coming from the last 
scattering surface, by the hot electrons population of the 
ICM. The photon energy variation, which is caused by the 
scattering process, can be expressed as CMB temperature 
variations 

ATsz = yTcAIBf{x){l + Sn{x, e,)) + ATfe^n (1) 

where 

y ^ J (^edTe = j {^~~^ arriedl oc j Pedl (2) 

represents the comptonization parameter, x = 
{hv) / [ksTcMB) the dimensionless frequency, h and 
ks are respectively the Planck and Boltzmann constants, 
TcMB, rUe and ctt, the CMB temperature at z = 0, the 
electron mass at rest and the Thomson cross section, 9e 
represents the dimensionless thermal energy of the ICM, 
Te is the electron optical depth. The parameters rie, Te, 
and Pe are the electron number density, temperature, 
and pressure of the ICM, Sn{x,9e) — fn{x)9^ / f{x) is 
the relativistic correction term that accounts for the 
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Table 1. Parameters of the electron number density profiles. 



Cluster 


rieo 








del 


Oc2 




/ 






(10-2cm-3) 








(arcsec) 


(arcsec) 








A1413 


3.89 ± 0.54 


0.535 


± 


0.016 


6.7 ± 1.4 


40.1 ± 4.1 


0.760 


± 


0.020 


A1689 


4.15 ± 0.31 


0.871 




0.040 


21.6 ± 1.0 


104 5 ± 5 3 


0.870 


± 


0.010 


A1835 


11.3 ±0.4 


0.802 


± 


0.015 


9.3 ±0.2 


63.8 ±1.6 


0.940 


± 


0.001 


A2204 


20.4 ± 1.1 


0.716 


± 


0.028 


7.5 ±0.3 


67.6 ±1.9 


0.959 


± 


0.004 


A2261 


4.07 ±0.59 


0.631 


± 


0.024 


10.2 ± 1.8 


39.1 ±5.9 


0.760 


± 


0.050 


MS1358.4+6245 


9.63 ±0.79 


0.676 


± 


0.017 


3.3 ±0.2 


37.0 ±1.8 


0.934 


± 


0.003 


RXJ1347.5-1145 


28.5 ± 1.4 


0.632 


± 


0.009 


4.0 ± 0.2 


23.3 ± 1.6 


0.942 


± 


0.004 


ZW3146 


16.9 ±0.3 


0.669 


± 


0.005 


4.4 ±0.1 


25.8 ±0.6 


0.882 


± 


0.004 



thermal energy of the electrons involved in the scattering 
processes, where f{x) = x[{e^ + l)/{e^ — 1) — 4] is 
a dimensionless quantity that describes the spectral 
signature of the effect, and the subscript n indicates the 
maximum order of the relativistic correction (n = 4 in 
this work. Nozawa et al. 2005,) . The last term of Eq. [l]is 
the kinematic component of the SZ effect, which contains 
the contribution from the bulk motion of the electron 
population with respect to the last scattering surface 
reference frame. For the purpose of this paper, this term 
is omitted, assuming that it is disentangled from the 
thermal component by multi-frequency observations, 
together with the signal from the primary CMB emission. 

The SZ effect is redshift independent and, for this rea- 
son, it is possible to detect distant clusters without any 
existing X-ray or optical observations. This is the case 
of the ongoing ground-based experiments such as SPT 
(|Ruhl et al. 2 004). ACT (Kosowsky et al. 2003D, and the 
all sky survey like Planck (jPlanck Collaboration 2011aP or 
the upgraded MIT O (|De Petris et al 2007^ and OLIMPO 
(jMasi et al. 2008^ with new spectroscopic capabilities and 
the proposed 30-m diameter C-CAT ( |Sebring et al. 2006D . 
However, some assumptions on cluster physics still have 
to be made in order to directly extract cluster observables. 

Estimates of cluster's total mass can be derived 
by SZ observations when X-ray or lensing mea- 
surements are available or by empirically calibrated 
scaling relations linking the SZ flux to the total 
mass (e.g. IVikhlinin et al. 20091 lArnaud et al. 20101 
IPlanck Collaboration 2"0TTbl IComis et al. 201ip . Total 
mass can also be determined by SZ observations 
alone when applying thermal energy constraints 
(jMroczkowski et al. 2011|) . 

To accurately reproduce the gas inside the cluster, an 
ICM universal model is mandatory (e.g. |Nagai et al. 2007 
lArnaud et al. 2010)) . In this paper we confirm that the 
simple isothermal 6eta-model is clearly an inappropriate 
cluster representation for total mass recovery by SZ obser- 
vations, particularly in the presence of relaxed cool core 
(CC) clusters. These objects show a well studied peaked 
density profile with a temperature decrement in the core 
region (j Jones fc Formarn984.) . In the local universe this 
class of clusters is observationally a significant percentage 
of the total cluster population (|Eckert et al. 20l"T|) . Even 
if the X-ray estimated CC fraction is biased by selection 



effects in flux-limited samples, recently a 35% of clusters 
have picked up in the SZ high signal-to-noise ratio Planck 
early cluster data-set (jPlanck Collaboration 2011c|) are 
CC clusters. Large scatter in mass estimates of CC clus- 
ters has been highlighted previously using numerical simu- 
lations by Hallman et al. (2006) and Hallman et al. (2007). 

We investigate the bias on the mass in a limited sample 
of eight nearby (0.1 < z < 0.5) and high-mass (M > 10"'^^ 
Mq) CC clusters observed by Chandra. The SZ maps 
of these clusters, which are expressed in thermodynamic 
temperature units and convolved with several instrumen- 
tal beams, are dealt with by applying the isothermal beta- 
model. The total mass is derived in three different ways: by 
assuming hydrostatic equilibrium and a fixed gas fraction 
and by applying a self similar scaling relation. To focus 
only on the consequences of the employed ICM model, in 
our analysis we neglect all the contaminants present in the 
sky by assuming in this way the best situation to recover 
cluster total mass. 

In Sect. [21 we discuss the electron number density ra- 
dial profile and follow self-similar studies to characterize 
a universal electron temperature radial profile of a lim- 
ited sample of eight CC clusters observed by Chandra. In 
Sect. [3] we generate maps of the SZ effect in thermody- 
namic temperature. In Sect. 2] we evaluate cluster total 
mass under different sets of assumptions. The bias on the 
recovered mass is described in Sect.[5l which discusses the 
main contributions. Conclusions are summarized in Sect. 

13 

2. Electron number density and temperature 
profiles 

A general parametric model for the cluster atmosphere 
must be defined to forecast the shape of cluster SZ signals 
in matched filter techniques for detecting clusters in blind 
surveys. ACT has detected new clusters assuming a two- 
dimensional Gaussian profile as filter (ISehgal et al. 20101 ), 
while SPT has detected a projected spherical beta profile 
(jVanderlinde et al. 2010p and Planck a universal pressure 
profile (jMelin et al. 201ip . 

The approach for determining the total mass cluster 
can be different. High-quality X-ray data allows an ac- 
curate modeling of cluster morphology, but in the case 
of low angular resolution and/or low signal-to- noise ra- 
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Table 2. CC galaxy clusters properties used in the analysis to generate a universal Tf. profile for this class of clusters. 



Cluster 


z 


Da 




6*500 




rj-i C 


name 




(Gpc) 


(kpc) 


(arcsec) 


(kcV) 


(keV) 


A1413 


0.142 


0.52 


1195 ± 232 


321 ± 62 


6.6 ± 0.6 


7.3 ± 0.2 


A1689 


0.183 


0.63 


1402 ± 260 


377 ± 70 


8.7 ±0.9 


10.0 ±0.3 


A1835 


0.252 


0.81 


1439 ± 414 


387 ± 111 


10.5 ±1.0 


8.4 ±0.2 


A2204 


0.152 


0.55 


1796 ± 320 


483 ± 86 


11.3 ± 1.8 


6.5 ±0.2 


A22Q\ 


0.224 


0.74 


1201 ± 168 


323 ± 45 


7.0 ±1.0 


7.2 ±0.4 


^^51358.4 + 6245 


0.327 


0.97 


1633 ± 885 


439 ± 238 


8.4 ± 1.1 


8.3 ±0.6 


/?XJ1347.5 - 1145 


0.451 


1.19 


1734 ± 170 


466 ± 46 


14.8 ± 1.5 


13.5 ±0.5 




0.291 


0.90 


1804 ± 344 


485 ± 92 


8.7 ±0.4 


6.6 ±0.1 



" IMorandi et al. 20071 
Temperature scale calculated by a weighted mean of Bonamente et al. (2006) data. 

Electron temperature of the cluster obtained by fitting the isothermal beta-model to X-ray data (jBonamente at al. 2006|l . 



tio a simple isothermal beta-model is still applied (e.g. 
[Marriage et al. 2010|[5ayers et al. 201ip . 

In this work we analyze this model 
(jCavaliere fc Fusco-Femiano 1978|) , which is based 
on the very general assumption that the electron temper- 
ature Tf. is constant along the whole considered cluster 
radial extension and that the electron number density 
follows a spherical distribution as 



(3) 



nejSo{r) = 1 ± 



where rigo is the central electron number density, Vc the 
core radius, and /3 the power law index. The subscript 
ISO indicates, hereafter, the isothermal beta-model. The 
proved inadequacy of this model is compensated for by 
the advantage of extracting a simple analytic expression 
for the y parameter along the off-axis angular separation. 



(4) 



yiso{Q) = yo 1 ± 



where 



2/0 = "-eO 



rUeC^ 



r(f/3) 



(5) 



and Oc = Tc/Da, with Da the angular diameter distance, 
which has been calculated for each cluster by using 



Da = 



Ho{l + z)Jo E{z'y 



dz' 



(6) 



where Hq is the Hubble constant and E{z) — 

[nMil + zf + ^a]^^'^- We adopt a kCDM cosmology 
with i/o = 70 km/s/Mpc, VLm ^ 0.3, and Vtp, = 0.7. 

It is easy to find clusters that are not relaxed or that 
display structures that are very difficult to describe with 
this model. Therefore, it is realistic to assume that many 
newly discovered clusters in blind SZ surveys exhibit such 
significant deviations as well. To explore a particular ICM 



gas morphology, we focus on CC clusters. We started an- 
alyzing a small sample of eight objects extracted from the 
Chandra dataset investigated in Bonamente et al. (2006). 

The study of a central region, commonly known as the 
core region, has been challenged by high-resolution numer- 
ical simulations ([Navarro et al. 1995 . Borgani et al. 2004 



[Kay et al. 20041 [Nagai et al."2007, Henning et al. 2009^ 
These works lead to an agreement on whether there is 
a cooling core in the very central denser gas region (r < 
O.lrsoo) of some clusters, as well as a slower decline in the 
temperature at large radii (r > 0.2r5oo). As usual we refer 
to rsoo as the radius of the cluster that defines a volume 
with mean density 500 times the critical density pcrit at 
cluster redshift. The choice of rsoo is motivated by simu- 
lation results from Evrard et al. (1996) showing the gas 
within this radius relaxed and in hydrostatic equilibrium. 

Moreover, many observational studies in X-rays have 
shown that the X-ray surface brightness, hence the under- 
lying density, cannot be represented correctly by a beta 
profile. A second component should be added or a peaked 
central part introduced in order to properly fit the ob- 
servation. The observed deprojected density profiles are 
peaked for CC systems and flatter for morphologically dis- 
turbed clusters. 

The rie cluster profile can be described by a double 
beta-model profile 



/ 1 + -^ 



i/3 



(1-/) 1 + 



' c2 



(7) 



where the parameters' data have been taken from 
Bonamente et al. (2006) and adapted to this work to have 
a symmetric standard deviation (D'Agostini, 2003). 

This distribution is a generalization of a double beta- 
model profile of the electron number density, developed 
by La Roquc et al. (2005), but instead using the same /? 
parameter for both the central region and the outskirts, 
as in Bonamente et al. (2006). The Vd — 6ci/Da and 
i'c2 = ()c2IDa are the core radii of the inner and outer 
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0.4 



Fig. 1. Radial electron temperature profiles of a 
Chandra-selected sample of CC galaxy clusters 
HBonamente at al. 2006\) and the best fit (solid line) 
of the electron temperature profile with the la error 
(dotted lines) proposed in this work. Temperatures and 
radii are expressed in terms of Tx and rsoo , respectively, 
which are used as scale quantities throughout this work. 



distributions, and / is a parameter defined between and 
1 that represents how the core region dominates the outer 
region. These parameters, together with Ueo, are taken 
from Bonamente et al. (2006) and summarized in Table [T] 
for the selected clusters. 

We describe the temperature profile as 



e,CC 



(r) 



T 



X 




for < 0.1 

''500 

for > 0.2 

rsoo 



(8) 



where A1.2 are determined by fixing the position at which 
the two power laws, described by the mi and rn2 param- 
eters, intersect each other, as explained in Appendix 1X1 

The search for a universal temperature profile in the 
cluster halo region has been a target of several works 
based on observations (e.g. IMarkevitch et al. 1998 
IDe Grandi ^ Molendi 2002[ [Zhang et al. 2004 

IVikhhnin et al. 20051 IVikhlinin et al. 2006 

[Sanderson et al. 2006J [Zhang et al. 2007 



[Pratt et al. 2007. ,Zhang et al. 2008 ) and hydrodynamical 
simulations (e.g.. ILoken et al. 2002i Borgani et al. 2004} 
[Kay et al. 2004[ IPiffaretti fc Valdarnini 200 81. 

The profile, proposed in this paper specifically for CC 
clusters, follows both the central drop and the outer de- 
cline of the gas temperature. The function is formalized in 
the log (r /rsoo) — log {Te{r)/Tx) plane on which the power 
laws of Eq. ID become linear functions (see Appendix IXI for 
a complete treatment). To fix the parameters Ai and A2 
of the radial electron temperature profile and to confirm 
the power laws indices mi and m2, we fit a co-adding 
of Chandra electron temperature normalized to Tx data, 
of the CC clusters selection with the proposed function. 



Tx represents the average temperature value in the range 
(0.1-;-1.0)r5oo, and it is used to scale each cluster, in order 
to fit the universal temperature function to the measured 
profiles. In Tabled we report the cluster redshift, z, and 
the angular diameter distance. Da- The scale radius rsoo 
and temperature Tx are also collected. The chosen radial 
range, which is used to calculate the scale temperature Tx, 
corresponds to a cut in the central region (r < O.lrsoo)- 
Obviously this value cannot be compared easily with re- 
sults of other works because it strictly depends on the 
data radial extension. In fact, an important source of bias 
is the temperature definition ([Vikhlinin et al. 2006). Here, 
we use the spectroscopic temperature Tx- In Figure [1] the 
temperature data of our cluster sample with the best fit 
are plotted. By following Appendix [21 the resulting Tf, 
profile parameters are Ai = 2.41 ±0.14, A2 = 0.55 ±0.10, 
mi = 0.38 ± 0.02 and ms = -0.29 ± 0.11, which uni- 
vocally define the Te(r) function. We note that, even if 
the power law that describes the outskirts of the cluster 
temperature distribution suffers larger uncertainties, mi 
and m2 are both compatible, within one standard devi- 
ation, with estimates available in the literature. For ex- 
ample, Zhang et al. (2008) find mi = 0.38 ± 0.04, for 
r < 0.2r5oo, in agreement with Sanderson et al. (2006), 
who proposes mi = 0.4, for r < O.lrsoo- For radii larger 
then 0.2r5oo, Zhang et al. (2008) fitted a selected sample 
of data from XMM-Newton finding structurally similar 
behavior to ours with m2 = —0.28 ± 0.19. 

3. Pipeline of cluster simulation 

The analysis reported in this paper can be summarized in 
the following steps: 

— construction of an SZ signal distortion profile 
ATszi^)), using the ICM information coming from ex- 
isting X-ray observations; 

— convolution of the cluster SZ map with several instru- 
mental beam profiles; 

— extraction of the ICM parameters as in the assump- 
tions of the isothermal beta-model; 

— estimation of the cluster total mass Mtot', 

— calculation of the bias on cluster total mass due to the 
incorrect description of the ICM. 

In the first step, we generate angular profiles of the 
SZ signal distortion ATszi^), assuming an observing fre- 
quency of 150 GHz. The electron number density n^. and 
temperature Tg profiles are constructed using the equa- 
tions presented in Sect. [5J which we assume to be a good 
representation of a cool core ICM. The angular profile of 
the comptonization parameter is then evaluated by pro- 
jecting the electron pressure profile on the plane orthogo- 
nal to the cluster line of sight. The SZ signal is obtained 
using Eq. [TJ 

The ATsz{0) profiles are then convolved by consider- 
ing three different instrumental beam profiles modeled as 
a Gaussian, corresponding to large single dishes (SPT or 
ACT) with FWHM = 1 arcmin, medium size telescopes 
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Table 3. Parameters of the single 6e to-model, estimated using the MCMC procedure described in the text, as best 
fit of the CC cluster maps. 



Cluster 


FWHM fov 




/3 






(arcmiri) 


{^IK) 




(arcsec) 


' " 1H:10 


1 


-273 ± 9 


0.731 ± 0.006 


66.1 ± 3.0 




4.5 


-224 ± 14 


0.743 ± 0.009 


81.7 ±7.3 




7 


-195 ± 11 


0.758 ±0.011 


96.4 ±8.0 


0/41 fiSQ 


1 


-554 ± 8 


0.979 ± 0.006 


80.2 ± 1.5 




4.5 


-430 ± 12 


0.994 ± 0.004 


93.7 ± 2.2 




7 


-592 ± 84 


0.909 ± 0.018 


61.4 ± 8.3 




1 


-709 ± 10 


0.951 ± 0.006 


56.8 ± 1.1 




4.5 


-468 ± 29 


0.980 ±0.013 


76.0 ±5.0 




7 


-818 ± 88 


0.870 ±0.010 


36.0 ±3.6 


( J I, A oori/i 


1 


-727 ± 8 


0.848 ± 0.003 


68.4 ± 0.9 




4.5 


-551 ± 18 


0.865 ± 0.007 


86.0 ±3.6 




7 


-514 ±39 


0.859 ±0.017 


86.9 ±9.0 




1 


-399 ± 10 


0.830 ± 0.006 


52.3 ± 1.7 




4.5 


-273 ± 17 


0.849 ± 0.008 


71.6 ±4.6 




7 


-323 ± 32 


0.816 ± 0.013 


55.5 ± 6.4 




1 


-304 ± 12 


0.857 ± 0.011 


49.6 ± 2.7 




4.5 


-245 ± 28 


0.856 ±0.013 


56.0 ±6.1 




7 


-1246 ± 116 


0.753 ± 0.009 


8.7 ± 1.2 


C?HXJ1347. 5-1145 


1 


-1683 ±15 


0.835 ± 0.002 


37.7 ±0.5 




4.5 


-893 ± 30 


0.860 ± 0.005 


62.7 ± 2.5 




7 


-739 ± 41 


0.873 ±0.011 


73.1 ±5.1 


Cizvi/3146 


1 


-645 ± 13 


0.839 ± 0.005 


42.4 ± 1.1 




4.5 


-404 ± 38 


0.857 ±0.013 


61.3 ±6.7 




7 


-585 ± 77 


0.801 ±0.011 


36.5 ±5.2 



(MITO or OLIMPO) with FWHM = 4.5 arcmin, and 
small apertures (Planck) with FWHM — 7 arcmin. 

The errors associated to the convolved ATsz{0) pro- 
files are treated as only due to instrumental noise. An 
optimistic choice of the sensitivity, for all the observato- 
ries, is 6 /.tK/beam, corresponding to the Planck channel 
at 143 GHz (jPlanck Collaboration 20TTc)) assuming the 
necessary integration time on source for the other experi- 
ments. Contaminants are not included in the study since 
we wish to assess our ability to extract the mass of the 
clusters under ideal conditions. 

To simulate the missing knowledge of X-ray observa- 
tional results, we ignore cluster morphology and assume 
the most general model for it: an isothermal &eto-model, 
that expressed in temperature is 

/ 5)2X2-2/3 

ATsz^ AT5zofl + ^] ■ (9) 

We apply the Metropolis Hastings (M-H) algorithm 
Monte Carlo Markov Chain (MCMC) procedure to fit this 
equation (after a convolution with the corresponding in- 
strumental beam) on the simulated profiles to extract the 
parameters ATszo, P a-nd Oc- For each cluster, at a fixed 
field of view {fov), we analyze the accepted set of param- 
eters derived by the MCMC procedure. The degeneracy 
among the extracted &eta-model parameters affects their 
uncertainties. 

We obtain the &eto-model parameter set that is most 
consistent with the ATsz{0) profile, given the assumed 



instrumental noise and beam sizes. All the parameters re- 
sulting from the MCMC analysis are collected in Table 
[3] for each cluster and for each fov. Figure [5] shows the 
electron temperature (top), number density (middle), and 
pressure (bottom) profiles for only the cluster ZW3146, 
as an example, of both the original CC ICM and the re- 
covered ISO model. The errors associated to the curves 
account for the la uncertainties on the parameters. 

4. Cluster total mass estimation 

We want to stress the consequences of the assumptions 
on the ICM physics when we miss X-ray information. A 
quantity, such as the total mass Mtot, can be biased by a 
different physical state of the ICM (i.e. mergers or cool- 
ing flow mechanisms). In particular we estimate the mass 
for both the ICM discussed templates (CC and ISO), by 
using the following different approaches: 

— hydrostatic equilibrium assumption for the cluster gas 
(hydrostatic equilibrium, HE); 

— gas fraction independence of cluster physical state 
(fixed gas fraction, FGF); 

— Mtot — y scaling relation (scaling law, SL), as derived 
in the standard self-similar collapse scenario. 

The masses are calculated, in particular, within a fixed 
integration radius rint (aperture radius), which we ar- 
bitrarily choose equal to the ^500 values as reported in 
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0.1000 



!^ 0.0100 ^ 



o 

„ 0.0010 



0.0001 



I 



0.100 



o 

> 0.010 



0.001 r 
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Fig. 2. Radial profiles of the electron temperature, 
(top), number density, (middle) and pressure, (bot- 
tom) for the cluster ZW3146. The red dashed curve de- 
scribes the CC template while the black solid curve repre- 
sents the ISO model (the shadowed regions define la un- 
certainties), with parameters extracted by MCMC analysis 
considering ATsz{0) profiles convolved with a beam of 7 
arcmin (FWHM). 



Morandi et al. (2007, see Table [2]). This choice is moti- 
vated by the need to fix an aperture radius within which 
to estimate integrated quantities. We point out that rint 
does not always correspond to the same overdensity, due 
to the different assumed ICM templates. It is clear that, 
hereafter, results associated to the clusters simulated in 
this work cannot be considered as describing the true ICM 
physics of the observed objects. We choose, however, to 
maintain the link with the "native" cluster in the name 
(NAME^CIname). 



4.1. Hydrostatic equilibrium 

The first approach, HE, assumes a spherical symmetry for 
the cluster, so that the hydrostatic equilibrium equation 
can be written (jSarazin 1988^ as 



dPgasjr) 
dr 



-pgas{r)G 



Mtot{< r) 



(10) 



where Mtot{< t) is the total cluster mass within ra- 
dius r and under the ideal gas assumption, Pgas = 
iPgaskBTgas)/itirnp), with Pgas = PempHe, p and pe are 
the total and electron mean molecular weights (i.e. the 
mean particle mass per electron in units of the proton 
mass rrip), G is the gravitational constant and Tgas = Ti = 
Te, the ion and electron temperatures respectively, be- 
cause the system is assumed to be in thermal equilibrium. 
To adapt this equation to the SZ physics, we substitute 
the gas pressure, which accounts for both electrons and 
ions, with the simple electron pressure P^, which causes 
the SZ effect, so we have 



The total mass can be derived as 

"91n(ne(r)) dlii{kTe{r)) 



91n(r) 



91n(r) 



(11) 



,(12) 



which reduces, in the simple case of the isothermal beta- 
model, to 



Mtot{< r) 



GprUp r^ -\- r^ 



(13) 



The previous equations yield Mtot,HE.cc and 
Mtot,HE,isOj respectively. If we also calculate the 
cluster gas mass by 



Mgas = PeTUp / n^dV, 



(14) 



it is easy to estimate the gas fractions, fgas,HE,cc, and 
f gas. HE, ISO- The plots in Fig. [3]refer, still as an example, 
to the cluster Glzwsue and show the radial profiles of the 
cumulative cluster total and gas masses for both the ICM 
templates and the gas fraction. The la lines are derived 
by considering the uncertainties of the ICM parameters. 

4.2. Fixed gas fraction 

The second approach, FGF, assumes that the gas fraction, 
estimated at rint , is independent of the cluster dynamical 
state, as deduced by simulations (e.g. IRasia et al. 20061 
ILau et aT 2 009) and as recently derived by XMM-Newton 
and Subaru observations ( [Zhang et al. 2010| . Under this 
hypothesis, we can derive the cluster total masses directly 
by gas masses as Mtot,FGF = Mgas /f gas- 
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Cluster 


ICM template 


Mtot,HE 


Mtot,FGF^ 




(FWHM fov) 


(10" Me) 


(10" Mq) 


A1413 


CC 




3.85 ±0.77 


8.73 ± 2.03 


CZai413 


ISO 


(V) 


6.59 ±0.60 


6.21 ±0.35 




ISO 


f4 5') 


6.67 ±0.59 


6.02 ±0.59 




ISO 


(T) 


6.66 ±0.58 


5.93 ±0.89 


A1689 


cc 




8.96 ± 1.97 


12.72 ± 2.33 


CZai689 


ISO 




13.56 ± 1.41 


9.88 ±3.12 




ISO 




13.61 ± 1.33 


9.14 ±4.08 




ISO 


(T) 


12.66 ± 1.33 


8.86 ±2.01 


A1835 


cc 




10.23 ± 2.30 


12.72 ± 1.18 


CIaiS35 


ISO 




16.21 ± 1.61 


10.08 ±0.33 




ISO 


f4 5'") 


16.32 ± 1.55 


9.13 ±0.94 




ISO 


(7') 


15.08 ± 1.40 


7.99 ± 1.44 


A2204 


cc 




12.91 ±3.30 


14.84 ± 2.96 


CIa2204, 


ISO 


{V\ 

\^ ) 


19.88 ± 3.27 


10.13 ±0.24 




ISO 




20.33 ± 3.05 


10.87 ±0.67 




ISO 


(T) 


19.89 ± 2.86 


10.32 ± 1.54 


A2261 


cc 




4.79 ± 1.11 


10.55 ± 3.48 


CIa2261 


ISO 




7.88 ± 1.09 


7.80 ±3.50 




ISO 




7.84 ± 1.18 


7.24 ± 1.22 




ISO 


(T) 


7.81 ± 1.08 


6.94 ± 1.09 


MS1358.4+6245 


cc 




8.09 ± 1.75 


10.23 ± 1.65 


C/jl/S1358. 4+6245 


ISO 




13.27 ± 1.77 


7.74 ± 0.63 




ISO 




13.17 ± 1.88 


7.27 ± 1.22 




ISO 


(7') 


11.90 ± 1.51 


6.12 ± 1.14 


RXJ 1347. 5- 1145 


cc 




14.66 ±3.02 


32.34 ±4.05 


CIrxJ1347. 5-1145 


ISO 


(1') 


24.45 ± 2.49 


25.04 ±0.48 




ISO 


(4.5') 


24.47 ± 2.54 


22.08 ± 1.30 




ISO 


(7') 


24.66 ± 2.49 


21.07 ± 1.84 


ZW3146 


cc 




9.06 ± 1.64 


17.88 ± 1.05 


^7/214^3146 


ISO 


(1') 


15.05 ±0.60 


13.77 ±0.58 




ISO 


(4.5') 


15.14 ±0.67 


12.72 ± 1.97 




ISO 


(7') 


14.44 ± 0.64 


12.03 ± 2.61 



" derived by Mgaa assuming fgas = 0.1 



Because the integrated quantities are calculated by re- 
ferring to a fixed radius Tint, defined in the previous sec- 
tion, it is evident that, for each cluster and for each con- 
sidered ICM model, we are dealing with different overden- 
sities. We calculate these overdensities by 



2G 



mm' 



-Mtotinnt)- 



The results obtained for the HE and FGF approaches 
are collected in TableSl In order to derive all the quantities 
corresponding to the ISO template from the parameters 
obtained with the MCMC procedure, we need to convert 
the ATszo values to central electron number densities. 
This is possible using equation 



AT. 



szo 



1 



(16) 



which is derived by expressing Eq.[5]in terms of ATsz- Of 
course this expression introduces a degeneracy between 
temperature and density, possibly producing the same 



ATszo value. For this reason the masses presented in this 
work have been obtained by fixing the electron tempera- 
ture to the Tx values (Table [2]) and assuming fgas =0.1. 



(15) 4.3. Mtot — y scaling relation 



The total cluster mass can be inferred by applying a self- 
similar relation that links it to the integrated comptoniza- 
tion parameter. The Mtot — Y scaling law is usually cali- 
brated for a fixed overdensity value, while in this work we 
are dealing with masses calculated at different overdensi- 
ties. For this reason in the scaling law, we make the de- 
pendence on overdensity explicit. Simple considerations, 
based on the assumption that cluster evolution is com- 
pletely determined by gravitational processes, lead to an 
easy scaling relation that connects the total mass of a clus- 
ter of galaxies, Mtot, to its temperature, Tg, considering 
an isothermal structure for the ICM. Following Kravtsov 
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Table 5. Total mass as derived by applying the Mtot — Y scaling law as in Eq. [511 



Cluster 


ICM template 


D\Y 


fgas 




Mtot,SL 


name 


(FWHM fov) 


(kpc^) 






(10" Mq) 


A1413 


CC 




100 ± 1 


0.23 ±0.05 


345 ± 69 


1.94 ±0.23 


CIauis 


ISO 




95 ±5 


0.10 ±0.01 


591 ± 53 


2.81 ±0.18 




ISO 


f4 5') 


92 ± 12 


0.09 ±0.01 


598 ± 53 


2.86 ±0.30 




ISO 


(T) 


90 ±9 


0.09 ±0.01 


597 ± 52 


2.81 ±0.26 


A1689 


CC 




177 ±2 


0.15 ±0.03 


477 ± 105 


3.34 ±0.38 


C/a1689 


ISO 




152 ±6 


0.07 ±0.01 


722 ± 75 


4.15 ±0.23 




ISO 




142 ±6 


0.07 ±0.01 


724 ± 71 


4.17 ±0.23 




ISO 


(T) 


139 ± 32 


0.07 ±0.02 


674 ± 71 


4.16 ±0.82 


A1835 


CC 




263 ±2 


0.13 ±0.03 


468 ± 105 


4.48 ± 0.36 




ISO 




186 ±7 


0.06 ±0.01 


741 ± 74 


5.05 ±0.25 




ISO 




171 ± 17 


0.06 ± 0.01 


746 ± 71 


5.12 ±0.51 




ISO 


(T) 


152 ± 28 


0.05 ±0.01 


689 ± 64 


5.00 ±0.85 


A2204 


CC 




271 ±2 


0.12 ±0.04 


337 ± 86 


5.28 ±0.72 


ClA220i 


ISO 




241 ±6 


0.06 ±0.01 


520 ± 86 


6.78 ±0.47 




ISO 


f4 5'") 


231 ± 15 


0.05 ±0.01 


531 ± 80 


6.89 ±0.59 




ISO 


(T) 


222 ± 31 


0.05 ±0.01 


520 ± 75 


6.88 ±0.99 


A2261 


CC 




134 ±2 


0.22 ±0.06 


388 ± 90 


2.31 ±0.38 


CIa22(,1 


ISO 




110 ±5 


0.10 ±0.02 


639 ± 89 


2.87 ±0.19 




ISO 


f4 5'1 


104 ±9 


0.09 ± 0.02 


635 ± 95 


2.93 ±0.29 




ISO 


(T) 


99± 16 


0.09 ± 0.02 


633 ± 88 


2.93 ±0.48 


MS1358.4+6245 


CC 




166 ±4 


0.13 ±0.03 


233 ± 50 


3.82 ±0.35 


CiA/S1358. 4+6245 


ISO 


\^ 1 


127 ± 11 


0.06 ± 0.01 


381 ±51 


4.67 ±0.45 




ISO 




113 ± 21 


0.06 ±0.01 


379 ± 54 


4.56 ±0.76 




ISO 


IT) 


96± 18 


0.05 ±0.01 


342 ± 43 


4.36 ± 0.76 


RXJ1347.5-1145 


CC 




1165 ±4 


0.23 ±0.03 


305 ± 63 


8.09 ±0.52 


C'iiJXJ1347. 5-1145 


ISO 


(1') 


719 ± 14 


0.10 ±0.01 


509 ± 52 


8.68 ± 0.40 




ISO 


(4.5') 


654 ± 41 


0.09 ±0.01 


509 ± 52 


8.86 ± 0.60 




ISO 


(7') 


649 ± 66 


0.09 ±0.01 


513 ± 52 


9.11 ±0.80 


ZW3146 


CC 




312 ±4 


0.20 ±0.04 


201 ± 37 


4.42 ±0.31 


Clzwi\i% 


ISO 


(1') 


224 ± 11 


0.09 ±0.01 


334 ± 13 


5.21 ±0.23 




ISO 


(4.5') 


208 ± 38 


0.08 ±0.01 


336 ± 15 


5.31 ±0.81 




ISO 


(7') 


205 ± 43 


0.08 ±0.02 


321 ± 14 


5.41 ±0.11 



et al. (2006) and Bryan & Norman (1998) and keeping the 
overdensity dependence explicit, we have 



1/3 



2/3 



To connect a spherically integrated quantity (e.g. clus- 
ter total mass) to one derived by a cylindrical integration 
(e.g. integrated SZ flux), we consider the spherical analo- 
gous to the SZ flux, Ys- This quantity is directly propor- 
tional to the cluster total mass as 



fcsO-T„ fgas Mi 



Ys = / UeTedV = 5-T„^,- 



gas-^'-l-tot 



(18) 



where Mtot is the cluster total mass inside a sphere with 
radius equal to r^nf, and Tmw is the gas mass- weighted 
mean temperature. 

By combining the previous equations, we get 

Mtot = 0.248[/3,,i?(z)2A,„i]"'/V|/'l0i4 Mq. (19) 



This means that, by estimating Yg, we can easily infer the 
corresponding total cluster mass. Furthermore, in order to 
derive the integrated SZ flux, we need to solve the equation 



(17) r = 2^ 



(20) 



where 9int = fint/DA and y{9) is extracted by the simu- 
lated SZ temperature decrement profiles (in the CC case) 
or calculated directly from Eq. |4j by considering the pa- 
rameters as derived from the MCMC analysis in Sect. [3] (in 
the ISO case). The two defined integrated Comptonization 
parameters {Ys and Y) are connected by the dimension- 
less quantity C = D\Y/Ys. 

We averaged the C factors, assuming different cylindric 
depths equal to 2, 5, and 10 Tint, obtaining C2 = 1.52 ± 
0.07, C5 = 1.85±0.14, and Cio = 2.00±0.19, respectively, 
converging to the value in Bonamente et al. (2008). 

The scaling law can be rewritten as 



A,;, 



'^'\D\Yf''>lQ^^ Mo, (21) 



which can directly give an estimation of cluster total mass 
Mtot,SL- All the quantities included in the scaling law are 
listed in Tabled 
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Table 6. Biases on the cluster total mass estimate for three different approaches as derived by three different fovs. 



Cluster 


FWHM fov 


MBhe 


MBfgf 


MBsL 


name 


(arcmin) 








C*/ai413 


1 


0.75 ±0.31 


-0.29 ±0.19 


0.58 ± 0.32 




4.5 


0.79 ±0.35 


-0.30 ±0.15 


0.59 ± 0.34 




7 


0.81 ±0.31 


-0.33 ±0.17 


0.61 ±0.31 




1 


0.56 ±0.31 


-0.24 ±0.14 


0.43 ± 0.24 




4.5 


0.59 ±0.33 


-0.30 ±0.15 


0.45 ± 0.28 




7 


0.50 ±0.30 


-0.35 ±0.20 


0.55 ± 0.45 


CI A1S35 


1 


0.63 ±0.31 


-0.26 ±0.07 


0.44 ± 0.16 




4.5 


0.64 ±0.28 


-0.32 ±0.09 


0.47 ±0.21 




7 


0.51 ±0.32 


-0.41 ±0.11 


0.45 ± 0.29 


CIa220A 


1 


0.59 ±0.32 


-0.25 ±0.15 


0.45 ±0.27 




4.5 


0.62 ±0.30 


-0.29 ±0.15 


0.50 ±0.27 




7 


0.63 ±0.30 


-0.31 ±0.18 


0.47 ± 0.33 


CI A22GI 


1 


0.75 ±0.30 


-0.23 ±0.25 


0.52 ± 0.36 




4.5 


0.74 ±0.33 


-0.28 ±0.29 


0.53 ± 0.40 




7 


0.69 ±0.32 


-0.30 ±0.25 


0.54 ± 0.40 


CImS1358. 4+6245 


1 


0.71 ±0.31 


-0.28 ±0.12 


0.50 ±0.21 




4.5 


0.69 ± 0.34 


-0.32 ±0.17 


0.55 ± 0.34 




7 


0.52 ±0.29 


-0.43 ±0.14 


0.47 ±0.32 


CIrxJ1347.5-1145 


1 


0.79 ±0.32 


-0.26 ±0.08 


0.51 ±0.17 




4.5 


0.77 ±0.30 


-0.35 ±0.08 


0.56 ±0.19 




7 


0.82 ±0.38 


-0.38 ±0.10 


0.61 ±0.26 


ClzW314e 


1 


0.66 ±0.32 


-0.28 ±0.06 


0.47 ±0.16 




4.5 


0.67 ±0.32 


-0.33 ±0.12 


0.53 ±0.27 




7 


0.62 ±0.33 


-0.40 ±0.14 


0.57 ±0.31 



5. Bias on the total mass 

The results clearly show that a general interpretation of 
cluster physics, as assumed in a pure self-similar scaling 
law like Mtot — Y, can infer a wrong estimate of cluster 
total mass. To emphasize and to quantify this point, we 
define a mass bias as 



MB = 



totJSO 



M, 



tot,CC 



My 



tot.CC 



(22) 



This quantifies the difference between the SZ derived 
mass, as results from the isothermal &e to-model assump- 
tion (ISO) and the X-ray derived mass (CC), which we 
deduce from the HE, FGF, and SL approaches. We no- 
tice no significant bias dependence on the beam size used 
to convolve the SZ signals, under the assumption of the 
presence of instrumental noise alone. Therefore, all plots 
refers to results obtained with a fov of 7 arcmin FWHM, 
considering Tx values as in Tabled] We calculate the mass 
bias for all the analyzed approaches. 

Since the l\Tszo parameter does not give us unique 
information on the physics of the ICM, we have to study 
different pairs of the parameters n^o and T^. Thus, we se- 
lect values of that describe a reliable range of electron 
temperatures (from 5 keV to 15 keV) and calculate the 
corresponding Tieo. The mass biases are plotted in Figure 
3^, where the three cases (HE, FGF, and SL) are shown 
all together in the plots corresponding to clusters CIaw?,^, 
CIa22M, and CI rxji^ht. 5-114.5, representing a wide span 
in the electron temperature values. Table |6] lists the mass 



biases, as estimated at the X-ray derived electron temper- 
ature {Tx in Tabled]). 

A check of the goodness of the procedure was done by 
simulating the observation of a cluster having an isother- 
mal 6eta- model ATsz{S) profile instead of a CC cluster. 
For this validation procedure, named Test, we used the 
same analytical expression of the SZ signal to extract the 
parameters. The assumed electron temperature is the one 
obtained by Bonamente et al. (2006) and reported in the 
last column in Tabled] In Figure HJd we represent the clus- 
ter total mass bias with these assumptions. While the bias 
is always zero in the HE approach, for the FGF and SL 
cases, as expected, we notice a mass bias dependence on 
Tg. Due to the degeneracy between the electron temper- 
ature and number density, the mass bias varies with rig. 
In all cases, it is worth noting that these biases nullify for 
electron temperatures equal to Teo values, thus proving 
that the method is not affected by systematics. 

The degeneracy between the ICM parameters, result- 
ing from yielding the same SZ signal, produces different 
trends on FGF and SL biases with electron temperature. 
For increasing values of electron temperature it underes- 
timates the Mgas- Considering the FGF approach, this 
implies an underestimation of the total mass, too. For the 
SL method this would instead produce an increasing trend 
in the mass bias with temperature because an Mgas un- 
derestimation corresponds to a decreasing gas fraction. 

The net result is a mass bias that is always different 
from zero for the HE approach, as well as for the FGF and 
SL ones, unless it is for an electron temperature suitable 
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Fig. 3. Radial profiles of the cumulative total (top) and 
gas (middle) cluster masses and the gas fraction (bottom) 
for the cluster Clzwsue- The plots refer to the two ICM 
templates (CC and ISO) discussed in the text, under the 
hydrostatic hypothesis. The ISO model parameters have 
been extracted by MCMC analysis considering ATszi^^) 
profiles convolved with a beam of 7 arcmin (FWHM). 



to correctly modeling the CC cluster as an isothermal one 
but generally different from the true ICM temperature 
(see the crossovers in mass biases in Figure S]) . 

The mass bias reflects a combination of different con- 
tributions (jPiffaretti fc Valdarnini 2008^ . In particular for 
this work, the bias can be reduced to an integration bias 
and a modeling bias, which are estimated and discussed 
separately. The integration bias refers to the mass bias 
due to the different integration ranges along the line of 
sight. An infinite integration range is considered for the 
analytic ATsz expression, while a finite range is assumed 
in the numerical projection method, which is limited by 
the missing knowledge of the electron pressure profile out- 



side the region defined by the data. The integration bias 
that can be evaluated by still applying the whole proce- 
dure to an ISO template (Figure HJ;) . The modeling bias 
reflects the dependence of the bias on the assumed model 
for the ICM. In order to highlight it and to cancel the 
integration bias, we apply the procedure to the CC clus- 
ter, recovered as an isothermal one. In this case both the 
cluster signals (CC mock data and ISO recovered data) 
are obtained by projecting the three-dimensional electron 
pressure profile. For disentangling the integration bias, the 
ISO profile is also integrated over a limited range along the 
line of sight. The HE bias ranges between 50% and 80%, 
always implying a mass overestimation that is indepen- 
dent of the electron temperature. 

6. Conclusions 

We studied the bias that affects the estimate of cluster 
total mass by SZ observations when an isothermal beta- 
model is assumed to describe the ICM physical properties, 
specifically in the case of CC clusters when X-ray and/or 
Icnsing information is missing. 

It is well known that, rather than the central 
Comptonization parameter yo, affected by the choice of 
cluster profile modeling, an integrated quantity, like the 
parameter Y, appears to be a more robust mass proxy. 
Nevertheless, we have shown that CC clusters can gener- 
ate observed y maps in which the ICM morphology could 
still be substantially hidden, even for the current most 
sensitive experiments operating from the largest available 
nim/submm telescopes. 

While a general assumption of cluster morphology is 
efficient at detecting them in blind SZ maps, the possible 
mismatch with the actual cluster profile results in a mass 
bias. In fact simple ICM models, like the isothermal beta- 
model, applied to SZ observations can wrongly estimate 
cluster total mass in the presence of peculiar ICM dy- 
namics as in CC clusters, which are studied in the current 
analysis, and mergers. 

We analyzed the mass bias as derived in a limited sam- 
ple of eight CC clusters observed by Chandra, both nearby 
(0.1 < z < 0.5) and with high mass (Af > lO" Mq). 
Under the assumption of an isothermal 6eto-model, the 
cluster total mass was derived applying three different ap- 
proaches: the hydrostatic equilibrium equation, a fixed gas 
fraction, and a self-similar Mtot — Y relation. 

Assuming we had no information from X-ray observa- 
tions, we reported the bias on the derived total mass as 
dependent on electron gas temperature. Only in the case 
of hydrostatic equilibrium does this bias appear almost 
constant for the considered clusters in the range of 50-80 
%. Incidentally, we notice that an electron temperature 
value exists for which the FGF and SL mass biases vanish. 
This could be the only case in which a simple isothermal 
6eto-model accurately reproduces the mass of CC clusters. 

The large biases on total cluster mass recovery in 
CC clusters represent another reason to definitely discard 
the isothermal 6eia-model for this purpose and to firmly 
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(a) (b) (c) (d) 




T, [keV] T, [keV] T, [keV] T, [keV] 

Fig. 4. Total mass biases estimated with different approaches, HE (black), FGF (red), SL (blue), for three representative 
clusters of our sample (from top to bottom: CIai689, CIa2204, o.'nd CI rxjisat. 5-1145 )■ The ISO parameters used for 
the bias estimation are associated to SZ signals convolved with a beam of 7 arcmin (FWHM). In panel (a) the total 
mass bias is shown, while in panel (b ) the reliability of bias estimation has been tested. The two bias contributions are 
disentangled: integration (c) and modeling (d) biases. Assuming no prior for T^, the mass biases are plotted versus a 
reliable range of electron temperatures (5-15 keV). 



support more sophisticated models, with universal pres- 
sure profiles fe.g. lArnaud et al. 2010|) . This is already em- 
ployed for modeling cluster atmospheres in almost all the 
present blind-survey data reduction (SPT and Planck), 
and it is planned in the next future for ACT observations. 
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Appendix A: Electron temperature profile 
function for cool core clusters 

Here we describe the simple approach used to obtain the 
electron temperature profile function applied in this work. 
The assumptions at the basis of this treatment are 



the knowledge of functional relations of the electron 
temperature in two different regions of the cluster (see 
Sect. [2) , and 

the approximate location of the maximum of the 
profile. 

Equation [S] can be also written as 



log 



7711 log ( — ^ 
^ ° \ 1-500 ; 

7772 log I 



qi for ^< 0.1 

(A.l) 
92 for > 0.2 

1-500 



which is clearly obtained by self-similarity studies of Tg 
profiles and by assuming a universal Tg profile describ- 
ing CC clusters. Thus, we need a function that follows 
those linear trends asymptotically. Following Eq. lA.ll it 
is simple to associate Ai,2 = 10^^'^. For the sake of sim- 
plicity, in what follows we call Y — log (Te(r)/Tx) and 
X = log (r/rsoo). 

Though several functions can satisfy these general 
conditions, we consider the hyperbole as the simplest 
candidate for yielding a continuum function. It is ex- 
tremely easy to find the hyperbole equation by know- 
ing its asymptotes, only if they are symmetric with re- 
spect to the coordinated axes. In order to be in this sim- 
ple case, we need to change our reference frame into a 
more convenient one, which results in translating the old 
reference frame into the intersections point of the two 
linear functions and by rotating it by an angle an = 
0.5 [arctan(777,i) -f arctan(7772)], which is the angle between 
the old X-axis and the bisectrix of the asymptotes (i.e. the 
new X-axis). The new asymptotes, which are now referred 
to the new reference frame, are 



TR 



7771 cosafl, — smafj 



cos aR + mi sm a n 
7772 cos aR — sin a n 



cos an + 7772 sm an 



X = niRiX = turX (A.2) 
X = mji2X = -777i^A(A.3) 



where the TR apex indicates that the equations are trans- 
lated (T) and rotated (R). Now that we have the two sym- 
metric asymptotes with respect to the coordinated axes, 
we can use the hyperbole equation to build up our func- 
tion, which reads 



Y 



TR 




(A.4) 



b = amn and c = VoM-fe^. Since 



where a ■ 

the hyperbole of Eq. IA.4I is now related to the wrong 
asymptotes, we need to inverse-translate it and inverse- 
rotate it. The resulting functions are 



Y 



TRR- 



-BBX ± y^BB^X^ - 4:AA{CCX^ - 62) 



2AA 



(A 



yTRR-^T-^ ^ -BB{X - Xp) 

2AA 
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± 



^BB^{X - Xp f - ^AA[CC{X - Xp f - b^] 
2AA 



+ Yp 



(A.6) 



where 



AA = 



cos^iap) - ( - I sin^iap) 




BB = 



CC = 




(A.7) 



and 



Xp = 



q2 - qi 



Yp = 



mi — 7712 

m\q2 — m^qi 



mi — m2 



The parameters mi, m2, qi, 92, characterize the 
asymptotes of the discussed hyperbole. The complete 
function can be obtained by substituting {X, Y) with 
(log (r /rsoo) , log (Te(r)/Tx)) and only considering the 
negative sign. We add a fifth degree of freedom identified 
with c to the analysis, which gives the focal point position 
of the hyperbole, and we arbitrarily consider in the anal- 
ysis c = 0.2, which is the one that gives a minimum %^ in 
the fit with the cluster electron temperature data. 



